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In Paper II [N. G. Phillips and B. L. Hu, this journal] we presented the details for the regularization 
of the noise kernel of a quantum scalar field in optical spacetimes by the modified point separation 
scheme, and a Gaussian approximation for the Green function. We worked out the regularized 
noise kernel for two examples: hot flat space and optical Schwarzschild metric. In this paper we 
consider noise kernels for a scalar field in the Schwarzschild black hole. Much of the work in the 
point separation approach is to determine how the divergent piece conformally transforms. For the 
Schwarzschild metric we find that the fluctuations of the stress tensor of the Hawking flux in the far 
field region checks with the analytic results given by Campos and Hu earlier [A. Campos and B. L. 
Hu, Phys. Rev. D 58 (1998) 125021; Int. J. Theor. Phys. 38 (1999) 1253]. We also verify Page's 
result [D. N. Page, Phys. Rev. D25, 1499 (1982)] for the stress tensor, which, though used often, 
still lacks a rigorous proof, as in his original work the direct use of the conformal transformation was 
circumvented. However, as in the optical case, we show that the Gaussian approximation applied 
to the Green function produces significant error in the noise kernel on the Schwarzschild horizon. 
As before we identify the failure as occurring at the fourth covariant derivative order. 

PACS numbers: 04.62.-|-v 



I. INTRODUCTION 



Many physically interesting metrics are conformally related to an optical metric. When the metric has an optical 
form, the Gaussian approximation Q for the Green function yields a closed analytic form. For a given noise kernel 
in an optical metric, one can take advantage of the simple conformal transformation property of the scalar field's 
Green function to compute the noise kernel for the corresponding conformally-optical metrics. In Paper II we have 
derived the regularized noise kernel in the optical metric by the modified point separation method and worked out 
two examples: hot flat space and optical Schwarzschild metrics. Hot flat space is related to thermal fields in the 
(spatially-flat) Robertson- Walker universe by a time dependent scale factor, and the optical Schwarzschild metric is of 
course conformally related to the Schwarzschild black hole metric. In this paper we present the details for computing 
the regularized noise kernel. 

According to the procedure outlined in Paper II, the main obstacle here is the subtraction of the Hadamard ansatz. 
The divergent Green function is defined in terms of the optical metric while the Hadamard ansatz in terms of the 
physical metric. We need to re-express the transformed optical metric in terms of the physical metric. The defining 
equations for the geometric objects on the optical metric are transformed to the physical metric and solved recursively. 

Now the Green function series expansion can be written solely in terms of the physical metric. While for the 
ultrastatic spacetimes considered in Paper II the Hadamard subtraction was straightforward, for spacetimes only 
conformally ultrastatic this subtraction is non-trivial. Before the subtraction can be carried out and the Green 
function regularized, we need to study the conformal transformation properties of the point separation objects used 
to define the Green function. As we will show, regularizing the divergent structure to sufficient order for the noise 
kernel is no small task, when we include the conformal transformation between the physical and the optical metrics, 
(e.g., the fourth order term in the expansion of the renormalized Green function has over 1100 terms). It is at this 
point that the symbolic computation environment takes over. 

Since all the series expansions used are recursively derived, as a check of the full expression it is sufficient that 
we get the correct results for the lower order terms. Indeed our general expression contains and confirms the Page 
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Q result for the vacuum expectation value of the stress tensor in the Schwarzschild black hole. With the general 
expansion of the renormalized Green function at hand, one can choose a metric and compute the coincident limit of 
the noise kernel as outlined in Paper II. In Section || we compute the noise kernel for spacetimes conformally related 



to ultrastatic spacetimes. In Section [II we specialize to the Schwarzschild black hole and discuss the significance of 
our results in Sec. IV. The Appendices contain formulas whose usefulness goes beyond the specific approximations 
adopted here. 

II. NOISE KERNEL IN CONFORMALLY-OPTICAL SPACETIMES 

For a static physical metric gab, we consider the conformally related optical metric 
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e'^'^gab (2.1) 



where the conformal factor e~^'^ is the space-dependent function such that grr — 1, *-e., the metric gat is for an 
ultrastatic spacetime. In general, we use the overbar to denote objects defined in terms of the optical metric gab and 
without for those in terms of the physical metric gab- For a conformally invariant field, the Green functions on the 
two spacetimes are related via 

G{x, y) = e''^(^)G(a;, y)e^"(^). (2.2) 
This is our starting point. In terms of the Gaussian approximation 

=, , . kA1/2 sinhKf 

GG^nss[X,y) ^ —r^- — ^ -, (2.3) 

OTT^r (cosh Kr — cos kt) 

where r = (2 (3)(t) ^ and K is the period of the imaginary time dimension, we have the expansion of the physical Green 
function as 



Gcauss = — Q- 2 ' S~~2 1 "H" + 757^ '"'^ + 2^'^ 

8(T7r^ Stt-^ 6 180 
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6CTAT2+4Ar'*)| +o(cri,AT5) (2.4) 



where to = u}{x) and oj' = uj{y). 

We proceed as in Paper II and subtract from Gcauss the Hadamard ansatz 

1 /2A^ 



S{x,v)^j^\^^— + awi+a^W2j+0{a') (2.5) 
to get the renormalized Green function 

In the Hadamard ansatz, the V{x, x') term is not included since there is no logcr divergences present in the expansion 
of the Gaussian approximation to the Green function. This can also be viewed as an extension of the Gaussian 
approximation. 



Now the situation here is different from the optical case as the divergent terms present (2.4) and (2.6) do not 
directly cancel. Much of the work for regularizing the Green function will entail showing that indeed the difference 
between these two divergent terms is finite and to develop this finite difference to sufficient order to compute the noise 
kernel. To this end, we write the renormalized Green function as 



Gron — Gdiv.rcn + Gfin — 7- — T^M^ (2-7) 



with 
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+ ^(-^-6aAr^+4Ar^)| (2.8b) 

W = awi+a'^W2 (2.8c) 

Both Gfin and W have weU behaved coincident hmits. As these functions stand, it is only the last one, W, that we 
can readily handle. Appendix E of Paper II gives the series expansion in terms of the world function a defined with 
respect to the physical metric gab- Since we want the noise kernel in the physical metric, it is the covariant derivative 
commensurate with this metric we must use when computing the noise kernel. On the other hand, the function Ggn 
is defined in terms of the world function and the VanVleck-Morette determinant of the optical metric (jab- Thus the 
coincident limit expressions for a and the series expansion for A 2 derived in the Appendices of Paper II cannot be 
used to determine the contribution to a series expansion of the renormalized Green function Grcn- 

In Appendix ^ here we show how to get around this problem. Both a and A 2 are defined in terms of covariant 
differential equations with respect to the optical metric. The conformal transformation properties of the covariant 
derivative are used to re-express these equations in terms of the covariant derivative commensurate with the physical 
metric. Then end point series solutions built from the physical metric world function a are found. The details, along 
with the found solution are collected in the Appendices of this paper. Using these results, we can determine the 
contribution from Gfin to the noise kernel. 

This leaves the first function defined above to deal with. The key to unlocking this term is to introduce the 
symmetric function 

S(x, y) = e-(-)a(x, y)e"(^) - a(x, y) (2.9) 
The important property of this function as shown in Appendix ^ is S '--^ ct^ as ti — s- 0. With this function, we have 

a (^Ai - A^j - S As 

Gdiv,ron = -Z — 2 7 (2-10) 

8 7r^ (7 ((T + 2j) 

To see that this is indeed finite, we use the expansions of 

S - (725.(4) (2.11a) 

A5 ^ 1 + (7A(2) (2.11b) 

A5 - l + crA(2) (2.11c) 

to obtain the leading order behavior 

1 -A(2) + A('^) + gW+^A(2)gW 

Gdiv.ren- g ^3 1 + a 5(4) ^ ' 

which is finite as ct — > and has the value 

[Gd™]-^(A(2)-A(2)-5W) (2.13) 



Using Eqns (Alia) and ( A16a), along with (C8) of Paper II, we get the explicit form of the leading order expansion 
tensors used above as 

.^(4) _ ^-.c^-dgab i^-cd9ab 'J^-p^'^ 9ab 9cd \ 

Kbcd - ^ +^ ^ (^•14a) 

Ai^,^ = ^ (2.14b) 
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t(2) _ ^;a^;fc , ^:afc ^■p^''^9ab , ^ip^Qab , Rgb -. , n 

^ab — a + a + 10 +10 (,^.i4Cj 



From these we can get the expansion scalars we need for (2.13): 

5(4) . 4pWp^I^r,l. = + -.i^^ + (2.15a) 

A(2) = 2pVA^'^ = (2.15b) 

A(2) = 2pVA^'^ = _ + ^ + ^■^v^-^iV'v'' + '^■,p<iV''v'' ^ P'^P'^Rpg (2.15c) 

3 6 3 3 6 
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where p° is a unit vector. From these expressions, we might expect there to be residual direction dependence for 
( ^l3| ). But when we substitute the expansion scalars into the cr — > value of Gdiv.rcm the direction dependences of 
the three expansion scalars cancel and we are left with 



[Gdiv,ren] = 77-^ (t^;/ " t^'^) (2-16) 
48 TT 



Using 



A5 



1 and [a] =0, we can immediately get the coincident limit of (2.8b) 



[Gfin] = j^^j—^ (2.17) 



48 e' 

Since [W] = 0, we readily obtain the coincident limit of the renormalized Green function 

which is the result derived by Page (Eq. (29) of [||). 

Now that we know we can regularize Gdiv.rcn: we turn to developing the series expansion of Gdiv,rcn to a sufficient 
order to compute the coincident limit of the noise kernel. We need the expansion 

1 ^-^(0) , , ^|r^(3) , ^2^(4) 



Gdiv.ren - (^div,ren + ^div,rea + ^div,ren + ^ ^div,rcn + ^div,rcn j (2-19) 

The computation of [Gdiv,ron], i-e. G^iil ,-cn' involves S to order cr^ and both and to order a. Thus to carry 

out the expansion (|2 . 1 9| ) we will need these functions expanded to order cr'* for S and cr^ for A^ and A^. These are 
done in Appendix \A\ For this section we will use expansions in terms of the scalar coefficients: 

S ^ ^(4) + ^1 5(5) _^ ^3 ^(6) ^ ^(7) ^ ^4 ^(8) (2.20a) 

A^ ^ l + aA(2)+ai A(3)+a2A(4)+aiA(5)+a3A(6) (2.20b) 

As ~ l + aA(2)+a5 A(3)+a2A(4)+at A(5)+a3A(6). (2.20c) 

(Recall the scalar expansion coefffcients are related to the tensor expansion cofficients via A^"^ ~ 2^pP'^ ■ ■ 
where in general we take -v/ctp" — c'"). 

Putting these in ( 2.1C| ), the expansion coefficients G^"^ ^.^^ are 

Gli^len = A(^) - A(^) - (2.21a) 

= A(3)-A(3)-5(5) (2.21b) 

Gfl^^.^^^ = A(5) - A(5) _ A(3) ^(4) _ A(2) ^(5) + (_A(2) + A(2) + ^t^' 



+5(4) |^_A(3) + A(-'') + 5(5)) - (2.21d) 
A(6) _ A(6) - A(*' 5(*' - A(3) 5(5) + 5(5) (-A(3) + A(3) + 5*5)) 
A(2) + A(2) + 5^) (5(^)' - 5(6))) - A(2) 5(6) 



g(**) = 

div,ren 

+5(4) (^_A(4) + A(4) + A(2) 5(4) + 5(6)) - 5(«) (2.21e) 

We have shown that for G^°], a direction dependence could arise, but when we use the values of the expansion 
tensors for E, A^ and A^, this direction dependence cancels. We expect this to happen because we see the coefficients 
that contribute to any given coefficient of Gdiv.ron are of a higher power than the G^".] coefficient in question. When 
going from the scalar expansion coefficient form to the tensor expansion form, the rank of the tensor is the same as 
the order of the coefficient. And each of the gI") is made up of higher order coefficients. But as we find by direct 
substitution from the expansion tensors listed in Appendix ^ for each case the direction dependence always cancel. 
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In Appendix ^ we give the complete form of the coefficients G^"^ i-^jj, in their corresponding tensorial form. This is 
one of our main results: The regularization of the leading order divergence of the Green function, when Green function 
has been computed in an optical metric conformal to the physical metric in which the problem is given. What is new 
here is that this regularization has been carried out to the order needed for the computation of the noise kernel. This 
analysis has been developed so as to take advantage of the symbolic computing potential of current workstations. 

The last remaining obstacle in computing the noise kernel comes from Ggn- As it stands, it is defined in terms of 
the optical metric, while it is the physical metric with which we nee d to ta ke the co varia nt deri vatives that determine 
the noise kernel. By using the function E and the series expansions (2.11b), (2.11c) and (2.11c), along with the series 
expansions 



-a(w+w') _ ,,,(0,a) 



+ \/aLU 



(ha) 



Gfin can be expanded in terms of the physical a as 



1 



Gfin 

O 7 

where these expansion coefficients are 



G 



(0) 
fin 



/^G 



(1) 
fin 



,(2,a) 



(2) 
fin 



(3,a) 



(4,a) 



'G 



(2.22) 
(2.23) 

(2.24) 



g: 



(0) 
fin 



'-'fin - 



,(04) 



g: 



(2) 
fin 



G 



(3) _ 



- f A(2) ,^(0^1) + c.(2a)) + ^ f 2 <5r(2) _ ^(o,2)\ 

6 V / 180 V J 

-fA(3)c.(0'i)+A(2)c.(i.i)+c.(3.i) 
6 V 



g: 



(4) 
fin 



K 

~6 



— ( 
180 V 

(aW c^(oa) + A(3) c.(i'i) + A(2) u;(24) + ^(4,1)^ 
180 V 

-2Jr(3) 



(1,1) 



A(2) 

,,,(2,1) 



-2,5t(2)^(0'1) 



,(0,2) 
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(4 ,5r(2)' c^("d) - 6 ,5r(2) c^(o.2) + ^(0,3)^ 



+ UJ 
(0,3) 



(2.25a) 
(2.25b) 
(2.25c) 

(2.25d) 



(2.25e) 



The explicit expressions for these expansion tensors are given in Appendix |C. 

With this we have regularized and expanded Gdiv.ren, and expanded in the physical metric Gfin, both to fourth 
order. The Hadamard ansatz function VK's series expansion is derived in Appendix E of Paper II, also to fourth 
order. From here we proceed as we did in Paper II for the optical metrics. To review, once a metric is selected, the 
component values of the expansions of Gdiv.rcn, Gfin and W are computed symbolically on the computer. Then Eqn 
(3.10) of Paper II is used to determine the component values of the coincident limit of up to four covariant derivatives 
of Grcn, along with the needed covariant derivatives of the coincident limits. From here the component values of the 
coincident limit of the noise kernel are computed. This procedure is adopted owing to the large size of the expansions. 
In fact, initial attempts to delay the specification of the metric until a full determination of the noise kernel ended up 
yielding a general tensorial expression of nearly 60,000 terms. By working out from the expansion tensors in terms of 
their component values, we found this had the added advantage of enabling one to study each of the separate terms 
that go into computing the noise kernel. 



III. SCHWARZSCHILD BLACK HOLE 



We can now turn our attention to a specific example. Consider a massless, conformally coupled scalar field on a 
Schwarzschild black hole with mass M. The line element is given in the usual coordinates — (r, 9, ip, t) 



ds' = 1 



1 - 



2M 



dr^ + 



sm 



(3.1) 
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where r is the imaginary time in the Euchdean sector of the spacetime. With the conformal factor 

e^- = 1 - — (3.2) 
r 

this metric is the physical metric corresponding to the optical metric considered in Paper II. As in that case, the 
imaginary time dimension has periodicity n — 1/4M, corresponding to a temperature associated with the Hartle- 
Hawking state. We use the scaled spatial coordinate x = 2M/r = Xj^nr. With this choice, the black hole horizon 
r = 2M is at a; = 1 while spatial infinity is at a; = 0. 

We now use the results for the expansion tensors above to compute the noise kernel coincident limit, along with 
the stress tensor itself. For the stress tensor: 

(Ta'') = (poo)diag{ (l + 2a; + 3a:2+4a;3 + 5a;^ + 6a:'^ + 15x^) /3, 
(l + 2.x + 3.x2) (l + 4.T3-3a;^) /3, 
(l + 2a; + 3a;2) (l + 4 a;''' - 3 a;^) /3, 

-1 - 2 a; - 3 a;2 - 4 a;3 - 5 a;-* - 6 + 33 a;^ } (3.3) 



where 



480 7r2 



(3.4) 



This result agrees with Page's Eq. (83) [||. In his work. Page showed the stress tensor must satisfy a functional- 
differential scale equation under conformal transformations. He then found, via trial and error, a general solution to 
this equation. This became the basis for his computation of the stress tensor in the black hole metric. 

In contrast, we have worked directly with the Green function and the point separation definition of the stress 
tensor. The conformal transformation properties of the geometric objects that go into the Green function are studied 
and the corresponding series expansions are developed. As can be seen from Appendices ^ and |^, these expansions 
can get quite involved even for just up to the second order terms needed for the stress tensor. Thus our agreement 
with Page on the stress tensor serves as an anchor for this work. Mo reov er, the methods used for deter min ing the 
series expansions are fundamentally recursive. The zero order result ( ^.Ig ) and the second order result (3^) are in 



agreement with known results. Since the terms needed for the noise kernel are generated recursively from these the 
symbolically computed results for the Schwarzschild noise kernel coincident limit should be accurate, up to the validity 
of the Gaussian approximation to the Green function. These results are 



2 

^ (219 + 876 a; + 2190 + 4380 a;^ + 7215 a;"* + 10464 a;^-|- 
756 ^ 

+ 16920 + 21424 a;^ - 2984943 a;* + 219140 + 197314 
+ 180260 a:" + 3292493 a;i2^ 

(137 + 548 a; + 1370 + 3860 x^ + 9005 x^ + 98432 x^ 

2268 

+225080 x^ + 408752 a;^ + 2553371 a;** + 1725900 x^ + 1206822 x"^ 
+761100 a;"- 4090521 a;^^) 



2 



2268 ^ 



37 + 548 a; + 1370 x^ + +2180 a;^ + 3965 x^ + 37952 a:^ 

-77240 a;*^ + 102992 a;^ - 11973349 a;* + 1817100 a;^ + 1721382 
4630140 a;" + 12756759 a;i2) 



(3.5a) 



(3.5b) 



(3.5c) 



o 2 

iV^^/ ^ -^-2^ (-219 - 876 a; -2190 a:^- 2140x^-495 a;'' + 2976 a;^ 
2268 ^ 

-10984 x^ - 49872 x^ + 1327551 x* - 230916 x^ - 50834 x^" 

+95356 x" - 774845 x^^) (3.5d) 

o 2 

iV^^g« = -^-2^ (-219 - 876 X- 2190x2 - 5500x3- 10575 x''- 17184 x'^ 
2268 ^ 

-19888 x*^ - 7200 x^ - 10917561 x* + 1056828 x^ + 999526 x^° 

+952012 x" + 11087083 x^^) (3.5e) 

o 2 

iV/e" = -^-2^ (41 + 164 X + 410x2 - 860x3- 4255 x"*- 50704 x^ 
2268 ^ 

-107048 x^ - 179440 x^ + 1023059 x® + 159708 x^ + 329206 x^" 
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+478972 + 1465003 x'^^) (3.5f) 
Knowing the component values, we determine the trace to be 

39000 n 2 

iV/,^ = '^'""Yr + ^) ("27 + 31 2:) (1 + x') (3.6) 

As with the optical Schwarzschild case, the trace under the Gaussian approximation fails to vanish, which it should, 
for the massless conformal coupling case we are considering. By taking the trace of the coincident limit expressions 
for the noise kernel above, term by term, we find that this arises from the non- vanishing of the fourth order derivative 
terms such as [G^cn;p^q'^]- This is what we had expected, based on our analysis of the optical Schwarzschild metric 
undertaken in Paper II. It is not from our implementation of the conformal transformation. 

Our result yields a zero trace at spatial infinity x 0, or r — s- oo. Using the fluctuation measure 



A a bed 

at r — > oo, we get 



{TabTcdj — {Tab) {Ted) 



{TabTcd) 



47V„ 



bed 



^Nabcd + {Tab) {Ted) 



abed : tttt rrrr 9990 rrrr tt99 rr99 

\ . J73_ 137 137 J73_ J73_ il_ 

i-^abcd ■ 2OO 200 136 136 104 



(3.7) 



(3.8) 



As we can expect for a Hartlc-Hawking state, this matches the thermal results of hot flat space we obtained in Paper 
II. 

We also report the magnitude of the error at the horizon: 

ahcd : tttt rrrr 9999 TTrr tt99 rr99 
g4 : 1904% 1904% 894% 18278% 1775% 1775% 

These results show the Gaussian approximation has completely broken down at the horizon. 

This is not to say we can draw no conclusions other than to discover the inadequacies of the Gaussian approximation. 
What is important is the finiteness at the horizon of both the noise kernel components and the error as expressed 
by the failure of the trace of the noise kernel to vanish. In our computation of both the stress tensor and the noise 
kernel, we have discovered finiteness at the horizon is a "fragile" property. By this we mean any small error in the 
symbolic code would result in a noise kernel that diverges as 2: — > 1. This has lead us to develop more than one way to 
determine the series for the optical metric geometric objects, just to test the symbolic code. We arrive at the results 
( |3.6[ ) using more than one computational path. In contrast if there were one single error in the code, the resulting 
noise kernel cannot be finite on the horizon. 

It is the noise kernel itself, via its trace, that provides a measure of the error of the Gaussian approximation, i.e., 
we self-consistently compute both the noise kernel and its error. Statements that address correcting the error of the 
Gaussian approximation also apply to the noise kernel. Correcting the Gaussian approximation will amount to finding 
the terms that need to be included such that it satisfies the field equation to fourth order in cr°. When this corrected 
approximation is then in turn used to compute the noise kernel, it will have to correct the current noise kernel results 
( |3.5f[ ) in such a manner as to exactly cancel the current trace ( |3.6| ). Hence the correction to the noise kernel will itself 
be finite at the horizon. With this in mind, we can conclude the fiuctuations of the stress tensor, as measured by the 
coincident limit of the noise kernel, are finite at the black-hole's horizon. This is one of the main lessons learned from 
the analysis of the noise kernel, as derived via point separation. 



IV. DISCUSSIONS 



Let us summarize our findings pertaining to two sets of issues: the range of validity of the Gaussian approximation 
and the results and usefulness of our program in spite of this approximation. 

Despite its success for the stress tensor calculation there is no compelling reason to expect the Gaussian approx- 
imated Green function to produce reasonable results for expressions involving higher order covariant derivatives, as 
in the noise kernel. Nonetheless the Gaussian approximation is very relevant because it contains the leading order 
divergence. This structure will remain even with a better approximation, while this leading order divergence must be 
regularized. This step is needed regardless of what form of the Green function one adopts - future improvements to 
the Gaussian approximation remains desirable, or, if the exact Green function is derived in the optical metric only. 
Our work lays down the structure and provides the details for its implementation. 
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Now for the successes and failures of the Gaussian approximation as appHed to the computation of the noise kernel. 
On the positive side our results for the fluctuations of the stress tensor of the Hawking flux in the far field region 
checks with the analytic results of Campos and Hu A fringe benefit is that we can verify our procedure by 
explicitly re-deriving the Page result js) for the stress tensor. We note that in Page's original work, the direct use of 
the conformal transformation was circumvented by "guessing" the solution to a functional differential equation. Our 
result is the first we know where the methodology of point separation was carried all the way through to the final 
result. That we get the known results is a check on our method and its correct implementation. On the negative side, 
our calculation shows that the fluctuations of the stress tensor based on the Gaussian approximation is unreliable 
in regions close to the event horizon. We show this by checking that the trace anomaly fails to vanish there. This 
result is not unexpected, as can be inferred from our findings in Paper II. Corrections to the Gaussian approximation 
need be introduced to improve the accuracy. One important result which may be under-appreciated is that we find 
a finite expression for the noise kernel. Even though the error is large, as long as it is finite, we know corrections to 
the approximations will themselves be finite. Hence, where the approximate noise kernel is finite, we can expect the 
full noise kernel to be finite. That the noise kernel on the horizon of a Schwarzschild black hole is finite is in itself 
a qualitatively significant result. This dispels claims to the contrary based on intuitive arguments or less rigorous 
calculations |, |, |, ||. 

Along the way to regularizing the Green function to fourth order, necessary for the coincident limit of the noise 
kernel, we have developed the series expansions of the various geometric objects that make up the Green function to a 
high order. Once the additional terms necessary to correct the Gaussian approximation are determined, it is a simple 
matter to use the work contained herein to compute these correction terms to a high enough order. In this sense, this 
work not only lays down the tracks and defines the steps, but also provides all the details necessary for implementing 
the point separation program to calculate the regularized noise kernel for quantum fields in curved spacetime. To 
end with a more practical note, for black hole fluctuations and backreaction calculations ^ which requires results 
from investigations like ours here, we can only reiterate what other researchers have said (lamented) before: A better 
approximation to the Green function is to be desired. 



Acknowledgments 

NGP thanks Professors R. Wald for correspondences and S. Christensen and L. Parker for the use of the MathTensor 
program. BLH thanks Professors Paul Anderson, Larry Ford and E. Verdaguer for discussions. This work is supported 
in part by NSF grant PHY98-00967. 



APPENDIX A: CONFORMAL TRANSFORMATION 



Since we compute the Green function in the optical metric conformally related to the physical metric we are 
interested in, we need to know how the geometric objects we use conformally transform. We denote objects in the 
optical metric with an overbar and the covariant derivative commensurate with the optical metric by a vertical stroke. 
The two metrics are related via 

gab = e~'^'^gab (AI) 

We develop the series expansions of the world function a, to eighth order, and A^/^ to sixth order. We also need the 
function 

S(x,?;) = e'^(^)+"(s')a-a, (A2) 

whose series expansion readily follows once the series for a is determined. The most straightforward method for 
determining these series expansions is to start by considering the conformal transformation of the differential equations 
each must satisfy. 

Considering first a, it satisfies, in terms of the optical metric, the equation 

a = \a\,a\P = ^5^^^|p^|,- (A3) 
Using (T|q = a a = d'-aj this equation transforms to 
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and we now have an equation purely in terms of the physical metric. We still have 

[a] = 0, [a,a] = 0. (A5) 

and we assume the series expansion 

+ ''pqrstuvw<^ <^ <^ (J <J a a a [Ab) 

To determine the expansion tensors, we proceed as we did in Paper II for A^/^ and substitute the expansion into the 
differential equation (A4) and collect terms by their order in a"". 
The order a term must satisfy 

2T(2),,-e2. |^Tp(2)^^^(2)^p^Y(2)^^Tp(2)p^^2T(2)„pT(2)f,) =0. (A7) 

This has the solution 

T(^)„,^^, (A8a) 

which can be seen via substitution. Since ct is a symmetric function, we use the results of Paper IPs Appendix B, 
which give the odd order expansion coefficients in terms of the even order one. Thus the third order coefficient is 

T^'\,c = = (A8b) 

With this, we have set up the recursion. Now it is just a matter to proceed as with the computation of the series 
expansion Van Vleck-Morette determinant as presented in Paper II. The results for the rest of the expansion tensors 
for a are: 

2LO'r(4) ■ 9cd ,n \ QabQcd ■„ /An \ 

e T^'abcd = (2w;qW;6 - W;q6) ^^uj.pU^'' (A8c) 

-^-^u;,,{io,,io-^P-u:,/) (ASd) 

+8uj.c^;def — ^-cdef) 
9ab 9cd „ ;p o ;P 

+9 UJ-p UJ-ef^ + 12 Ci^ip Re'^ f^) 
+ T^9abgcdgeSUJ.pUJ.q (w'^ W'"' - 3 W'^*) (ASc) 

+ 60W;ca;;de W;/g + 40W;ct^;d W;e/g " 20 LU -cd ; C f Q 

— 10 LO-c^jJ -de fg +'^;cdefg) 
-^lH^ {4:8uj,pUJ,eUJ: fU};gUJ-P ~ 27 UJ.pUJ.qRe'^fP.g 

— 66 UJ-pUJ-e ^^'^ -J g + 50W;p W;e'' W;/g 

— 156 UJ-pUJ-e^-J '^■,g^ + 72uj.e'^;pf ^-g^ + 22 LO-p UJ-ef UJ-g^ 
+9uj-pUj''^ UJ-efg + 7SiUJ-pUJ-e^;fg^ — 24 W;j,e "J^i/g^ 
-6W;pe W;/^g " 12^;^ W;e/g^ " 6 UJ -p U) ■ e f'' g 

6 W;p (^-e^ fg — 32 LO-p LO-qe Rf^g'' + 96uj-,p U>-q UJ-e Rf'^g^ 
-b2LulpLU,qeRf\^) 

9abgcdgef^:,p {12 UJ -q UJ -g LU'^ Uj''^ - lA UJ ■qUl''^ Ul -g^ 



8640 

-1 

+ 27 UJ,qgL0'P'i - 6 UJ,qUJ,gP'i + l%UJ,qUJ - «f 



10W;,W'Pw.g«+9W;/tJ;g« -36W;gW.gW'P« 
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+6a;;,a;'%-6a;;,a;;ri?/«'') (A8f) 



+ 2QuJ.cde^^;fgh - 60 tJ;,; t^;e/g/l + 30 a;;cd W;e/g?l 
+ 12 LO-c^-defgh — ^:cdefgh) 

(480 W;p U}.^e<-^;f <-^;g<^;h '^'^ 



120960 

— 1360 a;;p CO-e ^-J '^''^ ^■,gh + 308 Wjp U)'^ U>-ef ^■,gh 

+1936 Lo. p Lo. ^ ll). Lo.gfi — 200uj.p^uj.f^ uj.gh 

— 2080 LO-p LO-c UJ-J UJ-g UJ-h^ + 1552 CJ^e W;/ UJ-^pg UJ-^h^ 

+ 1048a);p W;e - 504 W;pe W; /g W; /j^ 

+384a;;)5W;e - 444w;p W;e^ a;;/^?^ 

+ 1584a;;pW;eW;/ W;g/j^ - 1128u;;e W;p/ 
-828W;p W;e/ t^:gft^ + 6 W;pe/ t^ig/i^ 

+ 135W;e/pW;g/,P + 48 W ;p W; g W ; / W ; / ft 

-264W;e'^;p/ W;//, + 72u).pOJ-efOJ-g''h 

+ 84:UJ.pef iO-g^h - 90tJ;e/pt^:g^/i - 36 LO-p UJ'^ U)-ef gh 
-bOAuj.pUJ.e(^-Jgh^ + 144W;pe W;/g?i^ 

-288a;;pa;;e W;/g^/i + 9QL0-peU).fgPh 

+264uJ.pUJ.eOJ.jPgh - 48uJ.pel^-/gh + 60tJ:pCJ.eyghP 
+ '.i6uJ.pUJ.efg^h + 12uJ.pUJ-ef''gh - 48 UJ-j, UJ-e'' f gh 
+ 36LO-pLJ-efg^h + '^■,p LO-e gh — 48 0.';,, '.^'-.c'' fgh 

-10Q8u;,pUJ.,qUJ,eRf%P-h - 912w;p Rge/.g uj-m" 

+ 180 U;:pU;.,gRe'^fP,gh- 336 U!;p U!;q Rre Rg'^h'' 
-8uJ.pRg''h'' {niiO-eOJ-qJ + 14:UJ-qef - 45W;e/q) 
+8 Rg'^h^ (232 OJ-^p (jJ-^qU-^eUJ-J — 271 CJjp CJje 

+82(i;;pett;;g/ - 120a;;pW;, W;e/ 
+56a;;pa;;ge/ + 15w;pW;e/g)) 

9ab 9cd gef /-, d n 

(192 a;;p W;/! 



241920 

— 96a;;p a;;g a;;g/( — 360LO-^pUj-qU!-gUj''^ u-h^ 

+464:UJ-pUI-qUI-g'' U)-^h^ — 4:08 U)-p W-q LO-g (jJ'P LJ-h'^ 
160u}-pUJ'P UJ-qgUJ-h'' + 316 UJ ■ p LO ■ g LO-qP to- 
+ - 64uj-pgUJ-qP UJ-h^ - 64uj-pqUJ-gP LO-h^ 

— b'lAuj.pljJ-qljJ-gU-hLO'P'^ + 9A8W-pU)-gU)-qhOJ'^'^ 

— 128 uj-pg uj-qh Lu'P'^ + 216a;;p w-qUj-gh oj'^"^ 

+ 168 UJ -^p UJ ; q UJ' UJ.^gh^ — lAAoj.pLU. qLU'''^ LU. gfi'^ 

— 300a>:p LO.^qP bJ-^gh' — 40a>:p LU ■ q LU' LU ■ g''' h 

+ 184LJ-pLJ-qOJ'P OJ-^g'^ h — ^288 (jj-p LO-q UJ-g OJ-^h^'^ 
-96 W;p UJ.^qg CJ-^h^'' + 672 W;p UJ.^q UJ.^g LO -^h^^ 

— 26 A uj.^puj.^q gU-h^P + 240a;;p a;;g w'^^/j 

-168 UJ-^pUJ-^qgLO'P'^h - ^2i^-,ptO-qLO-gh'"^ 
-132 CJ-p CJ-q W-gh''P + 66 W-p UJ-q UJ-g^h'^ 

+ MW;p!jJ-,(,UJ-,g''"'h + 6(jJ-,pUJ-.qUJ-^g^hP 

— 78 LO-p LO-^q UJ-g''''li — 48 lO-^p tO-q LO'^'* gh 
-3uJ,pL0-,qL0-,r (15i?g\«''' - 13RgP'ir,h 

+30 Rg\P-^ + 15 Rg^h"'" + 13 Rg''P''-,h) 

+ -80W:pW;,W;/ (5i?gV + 4 V^") 

+48 LO-,p uj-.q (3 uj.,rg Rh'"''' + 5 u:.,r uj.,g Rh^P^) 

32LO,pLO-,q {3L0.r LO'"" Rg^ +ALO.rg Rh'""')) 
9ab9cd9e}9gh , r. -p -o -r ai ;r [pq 

241920 (6w;rW'-^a;'^ W -67uj-rL0' lo'^^ 

+31iO-rLO'''LO'P'' + 102a;.^« W'f'' - eLO-rCO'^'"' + 27LO.rCO'P'''' 

+15a;;^a;='"'f) ' ' ' (A8g) 



where = denotes equality upon symmetrization. 



11 



By using the expansion 



e^"^ ^1 - V^y/ap^LO-p + ap^ p'' {uO-p UO-q + l^;pq) 
\/2 3 

+ — a^ pP p"^ p^ p'^ {uJ-pUJ-qUJ-r^;s +6 UJ-p LU-qLU-rs + il^;pq OJ-rs 

0-2 pP pt p^ p^ (tJ.p Wjg Wjr ti^;s W;t + lOu-^p LU.qLU.rl-^;st 



15^2 



+ 15 + 10 + 10 

+ 5 UJ-p LU-qrst +'^;pgrst) 



1 



yu 

+ 15 + 45 + 15 

^\pq ^-rs ^-jtu 

+ 20 LU -p LU -q LU -r ■,stu + GO LO -p UJ -qr iO -stu + W LU -pqr UJ -stu 

+ 15 LO-p LO-q LO-rstu + 15W;pqW;rst„ + 6 LO-p LO-qrstu ^■,pqrstu)) + O (^^^^ 



(A9) 



we can get the expansion of S. Multiplying together the above series and (A6) and subtracting cr, we determine the 
expansion 



„p „Q „r I v^(5) „p „q „r t . „p „q „r „s t u 

= Hqrs^^ ^ + ^pqrsf'^ <^ <^ cr + hpq^^t-^a^ a"" a a a a 



+^pqrstuv<^ <^ a <J a a + ^pqrstuvw<^ <^ a a a a a 



where the expansion tensors are 



(4) 
abed 

(5) 

abcde 
abcde f 



E 



5ah . , ^ Qabgcd -p 

3 6 

.9a6 /„ . \ , 9ab9cd p 
: (2 W;ct^;de + OJ-cde) + ^ ^;p^;e-^ 



3%/2 



3V2 



(AlO) 

(Alia) 
(AUb) 



.9afc 

30 



{u-ct^-.d^-et^-f + 2a;;c W;dt^;e/ + 7 LO -cd ;€ f 



+6 UJ-cl^idef +3 UJ-cdef) 
dab 9cd ,„ -r , rj -p 

+9 W;p tJ;e/ + 12 LU.p UJ.q i?e' /) 
+ l^&gabgcdgef UJ.pU.q {UJ'PUJ'''^ - 3 UJ''"') 



(Allc) 



.(7) 

'^abcdefg 



9ab 



{6uj-cOJ:dOJ-eUJ-Jg + 6 OJ -c OJ -.de OJ ■ f g + 3 UJ -c U) -d U! -e f g 



45%/2 

+ ^5^;cdt^;efg + 4 Ci^;c ti^;de/g + '^^■cdefg) 



gab gcd 

270^2 



(18 ljJ;p (jJ^e " 50 UJ ■ J' UJ ■ f g + 18 ; p W;e W;/ W;g^ 



-24a;;e t^^ip/ t^^ig^ + 68 a;;p UJ-ef l-^-.g^ + '^ll-^-^p^J'^ UJ-^efg 

—24 UJ-p uj-e^^-.fg^ + 24a;;pe uj.jgP + 6a;;pe 

+ 12W;pW;e// + 6uj.pU!;ef^g - & -p ■J' f g 
+ 27uj,pUJ,qRePf'^.,g- 12 UJ,p {2uj,qUJ,e-7uj,qe) RfPg'i) 
gabgcdgef (^8u},qUJ-P UJ.g" - 3uj,qPuJ.J - 9uj,qgUJ'P'' 



180 V2 



-4uj,q UJ.gP" - 2uj,q UJ''P\ +2uj,q UJ ,r Rg^"! 



(Alld) 



^abcdefgh = ^ iW;c W;<i t^;e W;/ + 3 W;c W;<i W;/ 

+69 UJ.cUJ-d^;ef ^-gh + '27 UJ-cd'^;ef'^;gh 
+2&U).(,^-d^-ei^;fgh + -c iO -de ■ f gh 

+AOuj-cdeUJ-Jgh + ISul-cUJ-dUJ-efgh + 53 Ul -cd UJ -e f gh 
+ 10 UJ-ci^;defgh + ^UJ-cdefgh} 

~^^^r (-"^^ {^^■,pei0.fgPh - 4:UJ.peUJ.jPgh + buj.pUJ-efgh^) 
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-Muj.pef {(^■,gh^ - ^■,g^h) 

+9 {lf)U!-efp^;gh^ + 16 t^;pe W;/gft'') 

+29cJ'''cJ;e/g/i - 28uJ-ei^-Jgh^ " '^0 LO -e UJ ■ f g^ h 

+ 16 LJ-eUJ-fPgh) - 24W;eW;p/ ( 19 tJ^g/i^ + 4 ) 

-8W;pe {25UJ-fP UJ-gh - 2lUJ-fgUJ-}J') 

+39 W;/ W^gft,'' + 12ijJ-j UJ;g^ h) + 266 LO -p LO'^ UJ ■ g f ^ : g ll 
+208 W;e UJ. f UJ.pg UJ-h^ + 18 UJ.pLU. eUJ. f to. g LU- h to'^ 
+Au}.pOJ.e^-f i^7uj'P UJ.gh - 16uj.gUJ.h^) 
+ 16 (11 LU-qUJ-e^-f — ^Auj-p LU-^^-gf +41 tjjpe W;^/ 

+3W;pa;;5 W^e/ + "^luJ-pUJ-qef + 30cJ;ptJ;e/g) ^g^/i'^ 
-312 UJ.p LO-q LO-e g'^ -h ~ 912 UJ.p Rqef^ -g UJ-h'' 

+ 180 LJ-pLJ-qRe^f'^-gh - 336 (jj-p UJ -q Ere Rg^ h^} 

^°i512o''^ {-256a;;pg - 300a;;pW;gP[j;g/i'' 

~2AUJ,pL0,qg (4C^;^f« + llC^;,,'?^ + 7 \) 

-2Auj,pUj,q {6uj,ghP'i - 3L^;g\« + U-J^ + 2^^^%,,) 

+8w;pa;'^ (20w;,g + 3uj.qUj.gh'' + 18 UJ.^qUJ.^g'^h) 

+24:UJ.pLJ.qLJ.g {2uj.,k.P'' + 3UJ-P''h) + 256uJ.pUJ.gLJ.qP UJ-k'' 
+ 16Lj]pLj]q {29uj.gPLJ,h'^ ~18uj,ghUj'^P'') 

+ 8 UJ.^pUJ.q {9UJ'P UJ''' UJ.^gh — 12 UJ.g OJ'P UJ.^h'' — 5 UJ.g uj.^h^'^'') 
+2AuJ.p UJ-qUJ-gUJ-h UJ'P UJ''' 

-3 uj.,p uj,q uj,r (15 RgPh"'"- + 45 RgPf^'" - 26 i?gP«^,,) 
+8w,pUj,q i9u;,rU;,gR,P'"' + 18 uj.,r g Rh^"'' + 16 LJ,rg Rh''^'^) 

+ 48W;pW;g ( 2 W ; ^ CJ ' ^ " 1 5 W ; / ) Rgh''} 

9ab 9cd 9ef 9 gh , j, -n -r c » or 
2520 ^iP^iqi^-.r^ — DW;rWW^ 

+ 17w./tj'«'' + 6tJ;^tj'f«'^) (Alle) 

Turning to the last expansion we need, we recall the Van Vleck-Morette determinant on the optical metric, A^/^, 
satisfies the differential equation 

Ai/2 (4 - Qct) - 2A1/2 ipfflP = (A12) 
Using the conformal transformation property 

□ct = e^'^ (Da - 2uj''Pa.,p) (A13) 

we determine A^/^ satisfies the equation 

Ai/2 (4 - e^" (Off - 2w;Pa;p)) - 2e2'^Ai/2 .p^^P = (A14) 
on the physical metric. Since we have the expansion for ct, we need only assume the expansion 

Ai/2 = 1 + A(2VV« + A(3),aV%'- + A(4),,a''aVa^ + A^l^aP a" a' 

+A(^;,,,„aVVVVV" (A15) 



and substitute this and ( [Aq ) into (A14) and solve for the expansion tensors. Following the now well defined method 
outlined above, they are found to be 

^ab = A'^ab + li^;a^;b + u;.,ab)~ j^i2u;.^pU;'P -u;.^pP) (A16a) 

^abc - A(fl~^{2uj..,aUJ,bc+UJ.,abc) + ^{^i^;p^-7-^;p''c) (A16b) 

^abcd = AI^I^+ ^{LJ..cUJ,d+u;.^cd) Rab+ ^{7uJ.^abUJ,cd + 6uJ;aUJ.,bcd) 
+ Y2Q i^;aUJ.bUJ.,cUJ.d + '^■^■ai^-b'^-cd + ^UJ;abcd) 
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+ 3 {Aoj-pP OJ-cd - 12 W;pc W;/ + SLO-p^cd) 

-12 oj.p {u-c 0J;d OJ'^ + 2 oj'P uj-cd - 2 oj-c W;/ - io-cd^ + 4 a;;/d) 

-2 (2u;;pW;d - W;pd) i?/ + 4 (a>':pu;:f, +a;;p5) -R/d'^} 

+ ^T7T^ {2 w.p w.q (6 w'f - 14 os^" + 3 i?f «) 
1440 

-14w;pa;'f + 5u;;pPw;g« + 4u;;p5 + Uoj.pLV-g'^P} (A16c) 

+ 15 W;ab W;cde + 4a;;Q LO-bcde + '2'tO-abcde) 
~T44 ^'^'"'■'^ '^;«^) (2 W;a W;6c + W;abc) -Rde) 

+2 (3a;;cdW;p^e + eWjp^Wjcde - SWjcWjp^de + 2W;j5^cde) 
+24 (Wjc U);pd W;e^ + 3 W;pc W;de^ - 5 LO-pc U>-d^e) 
+2ui.p {bbO-cdJ' +UJ.cd^e-'^'^^:/de) 

-24W;pW;c (w'''W;de + W;dW;e^) 

-Swjp (3a;;cdW;e^ + 3a;'''a;;cde - 10w;cW;de^ + 7a;;ctJ;/e) 

+4 (20 U),p UJ:gUJ.c + 22 LU,q LU.pc + LU,pqc + 10 LU,p Rgc) Rd^ e'' 

+4 i?//;e UJ,g + W^pJ + W^p (5 i?cd;e^ + i?c<i'^e " 6 i?/;de) 

-2 (2LU.cUJ-pd + 2 W;p W;cd + Moj.pcd - 15a;;cdp) -Re*' 

-5Rcd;e {2uj.,pLu''P - W;p^') - 5 (4W:p W;/ - a;;p^c) -Rde 
;p W;c + ljJ;pc) + 6 Rde'^ l^;pc} 

-^nlo^ {2a;;pa;;q (12a;'f a;;e« - 10a;;e^'' + 3iV'P%) 

— UtO-pUJ-q'^ U-e^ — 28u)-pU)-qeCO''^'' — 7 UO-pLO'^ iV-q'' e 
buJ.pP LO.q'^e + 6 CO-pe U) ■g'^^ + 4 Wjpg W'^'^e + 5 OJ;p OJ-^g'' e'' 

+W;pW;,«Pe + 3a;;p (w;qi?P«;e + 2w.,ei?^'^)} (A16d) 

AitLe/ = ^aLef + ^abcdef + 9ab^fdlf + QabQcd^-ff^ + 5a65cd5e/ A^^'^ (A16e) 



181440A^J^gy. = -22556 W;aW;6W;cW;dW;eW;/ - 67668 a;;aW;6W;ca;;dW;e/ 
-48444a;;aa;;bW;cdW;e/ - 10296a;;a a;;b W;cW;de/ 

-2612a;;Q6a;;cd'^;e/ - llQl&UJ-al^;bc'^\def - 108a;;aa;;6 W;cde/ 
+900 U;;a6c i^';rfe/ + 1332 Wjob t^;cde/ + 360 a;;6cde/ 
+ 180 W;abc(ie/ 

+90 (W;aW;6 + W;a6) i?cd;e/ " 5040 i?/e';/ -Rp6gc 

-90f?de;/ (56a);at^;6t^;c + 54LU.ai^:bc " t^;ahc) 

-30 (751 a;;at^;b'^;ct^;d + 1166u;;a u;;& W;cd + TStu-ab'^^cd 

+ lQ2oj.a<^;bcd — 3uj;abcd) Ref 
+ 105 {LL);ai^;b + UJ:ab) Red Ref 

-1260 (9 )RpegdRe''f" (A17a) 



120960A^J^ = 31672a;;pW;aW;6a;;cW;dW''' - 9116w;aa;;6W;cW;dW;p^ 
+41376 WipWjoWibW^Wjcd + 43936 a;;pa;;aW;6a;;cW;d^ 
-15568 + 4088 

+15152 

^■,a ^-fb ^-fpc + 32416 W;p ljj-a^;bc^;d^ 

-2664 LO 

;a'^;6'^;c'^;p^d + 4320 U^-p (^-a '^■,bcd 

+2256 u;;p a;;a tO;b ijJ;cd^ + 4896 u;;p to-a tO;b>jJ;c''d 

-1196 LO -pP UJ-ab^^;ed + 3056 CL);pa Cl);5c O);/ - 2952 UJ -a -be ^ ■,p'' d 

-2736t<;;o W;p^ W;6cd + 1152W;pW;a^ W;6cd + 3264 W ; « Wjpfc W; c/ 

+8256a;;pW;ab'^;c/ + 4800a;;aW;p6a;;c^d - h42Auj.pUj.ab'^;J' d 

-12 (jj. a '^■,b(^;p^cd - 'i08uj.pUJ'P UJ.abed + 1848 W.p W.^ W;6c/ 

-264 LO-pLO- a t^-M^d - 624a;;p a;;o W;6^cd + 1560 W;pa6W;cd^ 
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+ 1260W;a6p W;c/ - 3180W;pabi:^;/d - 156 W;o6 W;p^cd 
+ 10S U.pP ijj. abed + 1728W;paW;bc/ ^ 192 tJ^pa W;6/d 

-2112a;;po W;bPcd - 144 W;a Wjp^bcd + 216a);pa;;a6c/ 
+96w;p W;afc/d - 24a;;pa;;o5^cd - 528t<;;p Wjo^bcd 
+60a;;pPa6c<i 

+ 12 W;p (9 Rab;ccf + 4 Rab-c'd — Rab'^ cd — 12 RJ' -bed) 

-6 (10a;;pW'^ - 21a;;p^) Rab;cd + 12 (3W;pW;o + 16W;pa) Rbc;d^ 

-12 (5w;p a;;a - 6a;:po) -Rbc'^d - 72 (w;pa;;a + 3w:pa) Rb^-cd 

+3360cc';p LU.^a t^'^ -R&c:d + 288lJ:p UJ;a t^:5 -Rcd'^ 

+696 W-p LO-a W;6 -R/;d - 96 a;;a Rcd'^ OJ-pt, - 96 LO-a R/;d l^-.pb 

+3960J-p Rcd'^ U)-ab + SSAoj.pRJ'-dUJ-ab + 600 LO-pRbc 

-480 Rcd'^ UJ,pab - 1320 RP,d i^;pab + 210 i?6c;d W^p^ „ 

+630 i?ed'^W;abp + 1260 W;abp + 84W;p {R.d'" " ^ab 

-840u;;p (i?ed'« + 2i?c«;d) iigab^- 96 (a;;p W;, + u;;pg) i?aV;cd 

-12i?b2'c'';d {d:Ou).pUJ.,qUJ.,a + 222u;;a W;p5 - 112cj;qa;;pa 

+ 15W;pqa - 70W;, i?pa) + 348W;p (i?b«/;d + i?b"c^d) i?ga/ 

+24c^;p (11 J?eV=^ + 15 i?c«P^d + 15 i?,arb 

+2 (8300 UJ-p Ul-a^;b W'^ + 42 LJ-a W;b W;p^ + 1688 a;;p U)'^ W-ab 
+840J.p'' LJ-ab + 6504u}.pU}.a>jJ;b'' + 132 LJ -pa OJ-b'' 
-A2uJ.aL0-p^b + 1188 CJ;pCJ;ab^ - IQSuJ-pLO-Jb + ^3 ■p' ab) Red 

+24 (246 W;pW;aW;(,W;c + 19W;(jW;(, Wjpg + 304 

+11 W;pa W;bc + 44ct;;a W;pbc + 26w;pa;;abc 

— 49 W;a W;f,cp — 19 I0;pabc — l^;abpc + 28a;;abcp) -Rd^ 

+56 (37a;;pa;;g + a;;pg) RabR/d'' 

-1120 (a;;pa;;„-2a;;pa) RgbRc^d^ 

—4 {1844 W;p uj-qUj-a^-b + 2308 a;;o Ci^;b W;pq + 256 

-llAA U-pbLO-qa - 1208 CO ■ pU). a l^^qb + ^6 iO;pa '^■,qb 

-1520 + 316 

t^ipg '^;ob + 660 UJ-a^;pqb 
--lQ8uJ.qL0.pab - 18 UJ -^p LO .^qab - 420 Wjp Wjobg + 24a;;p,ab 

-7 (30 LO-r t^''^ + t^ir'^) -Rpoqb 
-12 (59cj;rtJ:a + 2W;ra) ^^pbg"" 
-4 (191 W;r W;a + 20W;ra) Rp'' qb 

- ((19a;;pa;;^ + 16a;;p^) RqabDRJ'd'' (A17b) 



120960A^^j,'^ = -18232 a;;p a;;a W;b w'^ w''' + 10364 W;pW;aU;;bt^'''t^:g'^ 

-8w;pW;g (520 w w ' ^ w ; a b + 35 1 8 w ; a ' ^ w ; b'' + 1353u;;aW;bW'J'«) 
+14u;;aa;;b W;p^ + 2720u;;pa;'^a;;q' W;ob 
+9888a;;p a;;a a;;q^ a;;b^ — 2632a;;pa;'^a;;ga'^;b' 

— 4:576 Ul.pU.qU-a^ i-^-b'' — 2368 U>. a ^^■,bl^;pq<^'^'' 

— 9152 LO-p Ul.^a i^;qb t^'^'' — 3AA0 CJ-p LV-^g U!-ab 1^'^'' 
+2700a);p W;a W'^t^;q'^b + 240oJ.pUJ.atO;btO;q'^^ 

— 3024 ct;;p a;;<j w'^ W;ab^ + 1104 W;p W;g w'^ (x;;o'b 

— 96 LO.pLU.qLU. a '-^;b^'' ^ 2712 LO -p LU .^q LO -a '-U'^'^b 

+98 w|pP W;af, + 408 W;pa W ; g W ; b^ - 2096 tO.pq LO.^a^ U-b'^ 

+784:LO.paOJ;qbi^'^'' " ^24:UJ.pq UJ-ab i^'^" " 84 a;;^ W;pf W-.^'b 

+756W;pW;a^a;;,«f,- 120CJ;„CJ;p,,W;q«P+444W;pW;abW;,'''' 

+1656w;pa;;g''a;;ab^ — 2316a;;pW;q^ W;ab'' — 648a;;pa;;q'^ W;a^b 

+504 U.p U);gP W-a'^b + 864 U-p U-ga Ul-b^'^ - 2928 U-p UJ-qa Ul-b''^ 

— 1944 LO-a OJ.pq LO'^h - 1104 U).p UJ-qa UJ'^'^ b + 30 UJ-p U'^ U).q\b 

-72W;pW;a W;g''b^ - 24W;pa;;a W;,''^^ + 108W;p W;ab^'' 

— 648uJ.^pUJ.^qUJ.^a^b'^ — 120 U) -^p UJ -^q UJ ■^a'^'^ b + 132 UJ .^p UJ .^q UJ'^'^ ab 

+105u;]pPaCO,g\ - 480a);paba;;,«f + 630 u;;pP« W;abg - 90a;;p,„w'P% 

+ 126u,pPu.g''ab + 192u.paUJ.g\f + 72W;pa - 96u.pgUJ'P''ab 

+108cj;pCj;g\bP + 48a;;pW;,\P6 - 12 W;p a;;,«Pab 

+ 18 UJ,p U.g (3 Rab-"" - 6 RaP-b" + 2 RaP'\ + 5 RP'^-ab) 

-12u,puj.g (w'P (26i?«b=' + 53J?„9;6) + 2a;;„ {2 Rb^'" + R^" .b)) 
-2Au.,p (RbP'" -2Rb'>'P - 11 i?f«;b) u;;ga 

+6u.p {{17 Rab'" - 18 i?a^b) t^;,^ + 14 {Rab'" - RJ;b) ^.q") 
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+2Au,pU;,qUJ,r (lli^cA'?''^ - 30i?a''«'-;6) 

+24a;.p (11 i?„V'P + 30 t^;,, + 42 w.p w., 

- (1580a;;pa;;ga;'Pw'« + 98cj.p w-^ - 35w.pPw.q« 
+964w;pu;;gw=f« - 28cj;pq U;'P« - '84cj;pa;;,«P) i?a6 

—4 (1470 Wjp Wjq (X';a + 101 CJ;g CJ'' Ct;;pa + 16a;;<j W;a W;p' 

+ 14a;;p W;o W;,' - 7W:pQCJ;g* + 736W;p W:a' " 2W;pga;;a^ 

+180w;,a;;pa'' + 24w;5a;;p«a - 210 cj;, iJfcf 

-4 (89W;p W;qti;;at^;fc + 28W;p a;;a - Muj-pa i^;qb 

+ 63 CJ;p UJ-^ab - 36 tl^ ; p a;;5a6) -^^ ' + 672 W;p W;g Rr a Rb^'''' 

— 8 (179 W;p W;g CO -r W ' " 165 Wjr LO'^ LO-^pq + 378a;;q W;p'" 

— ^ U) -^p UJ -^r LO -^q^ -\- ^ LO -^pr U) -^q^ — 175 LO-^pLO-^qLO-^f^ — 7 OJ -^pq CO -^f^ 

— 33uJ-ri^;pq^ — 105a;;pU;;r'^9 + 4:5ijJ.r<jJ;s Rp^q^) Ra'^b'' 
4 (1124a;;g W-^a l^;pr + 400 a;;p U>-^q UJ-^ra — 104a;;pq Ul-^ra 

+ 174:U);qCO;pr'a - 87 OJ-qUJ-s Rpar' - 87 UJ-qW-s Rp\a) Rb^'"' (A17c) 



362880A(^') = 4792 w;p o);^ tj'P - 3732 W;pW;gW'Pa;'«a;;/ 
+8808a;;pW;5 - iAoj-pU'^ ui-q'^ lo-/ 

-2568 oj.p tj.q W;/ + 1440 oj.p w'^ uj.qr 

+ 1248 CJ;p U,q U;/ UJ''''' - 828 IV.p UO,q IV'P uu,/'' 

+360 cj;p a;;r w'P'^'' + 35 W;pP a;;^'' a;;r'' + 84w;pga;;r'' w'^' 

+ - 64 U.pq W./ + 252 UJ.p UJ-q'^ UJ./P + 306 'uj-p Ul-qP W;/' 

+792CJ;pW;q/cj'«'^P + 162tJ;p'tJ;q 

+18a;;p (56<x';5 w''' + ISwjr Wig*" + 7 uj-^qUj.^/) -R^^ 

+108a;;pa;;g (a;;/i?f«'''-10w;^s'i?P''«") ' ' (A17d) 



APPENDIX B: EXPANSION TENSORS FOR Gdiv ren 
We present the explicit expressions for the expansion tensors for Gdiv,ren- The expansion scalars in the body are 

related to these tensors via G^"],,en = 2^P"' • • •P"''<^d?v,ren,al...a„• 
^?^il,^en = ^(^;p^-a;;pa;=^') (Bla) 

GiLc„,a = ^(2c.,pa;;,P-o.;p^) (Bib) 

^dii,ren,ob = (4 W;p W;a W'^ + 15 W;p - 15 W;p i?^^;;, - 4 Wja W;b Wjp^' 

+2w-p'' (jO-ab + 8w-pU}-a<JJ;b^ — 20 LJ^patO^b^ + 6 UJ-a i^;pb^ 

-12 LU. a LU.p'''b + 36LU.pLU.ab^ - 54(^;pCJ;a^6 - 12u)-^pab^ 
+3W;pa^6 + 9uJ:p''ab + 9W;a&p^ - 5 UJ-p LV'^ Rab 

+5 u),pP Rab- 10 u;,pUJ,^ RfP + 5 oj,pa Rb"" + 34 cJ.p u),q Ra^b" 

+ 10 LJ,pqRa\n/3Q0 

—gab (2 OJ;p UJ-^qLO'^ U)''^ + ^UJ-^p OJ'^ OJ-q'^ — 5 UJ-p^ UJ-q'^ + 16 W;p UJ-^q W'^' 

-Aw.pqU)'P'''-12w.pW.q''P'-60J.pW.qRP1)'/726 ' ' (Blc) 

<^div,ren,afec = (25 W;p W'^' i?afe;c + 50 W;p W;a i?6c'^ - 10 W;p W;a -Rfe^ic 

-IZOuO-pUO-qRa'^b^-c - 25 Rab;c(^;p'' - 40 Wjp W'^ Wjfcc 

-t-Ztfl/,) _ r,t P I,! I _ — Fi9 I,! I,! P I,! I ^ 



+ 40 W;aW;pPW;bc " 52 W ;p W ; W ; fcc 

-40a;;p W;aW;6W;c^ + 80 i?pa;bW;/ - 100i?a6;pW;/ - 40u)-a'^;pb'^;c 

+ 12oj.pU!-ab'^;/ - 40 iJpaqbjcW'^'' - 21 UJ-p Rab;c^ - 39u)-p Rab'^c 

+60W;p i?a^;6c + 20 LU -a tO -b l^-pc^ - ^0 U! -ab ■pj' + QO UJ -ab ■p'' c 

-10 UJ ■ p'' UJ;abc + 20a;;p U;;a tU-b/ - 1^0 UJ- pa Ul-bc^ - QOuj-pLO-at^-b^c 

+ 2S0 pa -b^ c - 80uj;aUJ:pb/ + 10W;at^;pfc^c + 50W;a W;p^hc 

-110 Ul-pUl-ab/ - 30W;p Wjab^c + 180 U)-pU)-a^ be + 50W;a Wjfecp^ 

+20a;;pa6c^ + 20uJ.pab^c - 10 a;;pa^6c - 20 UJ-p^abc 

+ 10 U)-abp/ - 20 LO-abp^c - 20 UJ-abcp^ + 50 UJ-p LO-a^ Rbc 

-25 U-pPa Rbc + 10 W;p Ra^ Rbc - 20 UJ.p W;a W;b -R/ + 70 W;a W;p6 i?/ 



16 



+ 60 a;;p OJ-ab R/ - 50 UJ-.pab i?/ + 20 UJ-abp -R/ - 10 OJ-p Rab Re" 
+20 iV,p LO,q UJ.,a Rb"/ + 20 OJ.,pq 

-480 a;;p uj,qa Rb'^J' + 20 uj,pqa Rb^c" 

— 80 CO-^paq Rb'^c' — 32 U!-p Rga Rb'^ c" 

— 32a;;p Rqra" Rb^ + 24W;p Rqar" Rb^ 

+32 a;;p - 16 c^;p Rqarb Rc^H /3600 

-5ab (15aj;pa;;q i?/''' - 100a;;pW;g W'"^ W;/ 
— 20 a;;p a;:c^ + 80 u ■ p u ■ q uj' ^ lo-c'^ 

+ 132 UJ.p UJ. qP UJ,c'' - 212 UO.pUO.qcUo''"'' 

— lOWjp W^Wjgc'^ + 10W;p W;c W;g^^ 

+25 Wjp^' Wjg^'c - 10 UJ.p + 30 UJ.pc - 70 tJ.p UJ.q UJ.J"' 

+10 uj.pq + 30 W;p w.q J'P''c + 10 a;;p9 w''"^c - 20 w|p W;,/'' 

+ 5U.pL0.qc'^P + lOW^pWjgV - 5W;pW;/«c + 25W;pW;g«/' 

+ 5 a;;p a;;g«Pc + 10 CJ;p a;;/q« + 10 W;p W;, LO'P i?c* - 5 W;p LO.qP Re" 

-W'cO.^p'oJ.qOJ.cR"'' +35u.,pU}.qcR"") /3600 ' ' (Bid) 



^div,ren,ofccd = {~960 U!;p UJ;a l^;b iO;c <^;d C^'" + 2520 UJ -p UJ -a OJ;b Rc" -.d 

- 4520 U ;p W ■ qW- a Rb"c' id + 960 W;a t^;fc W;c W;pP 

-1440W;pW;a W;6W'^a;;cd " 6300 W;p Rab'^OJ-cd 

+480 + 2760 

-2640 LU-pP UJ.abi^;cd - 48544 UJ -p UJ- a UJ-b'' U)-cd 
+7384 LO -pa LO -b^ LO-cd " 960 Wjp W;6 W;c 

+5040 a;;a iip6;c iv-d" - 6300 a;;p i?a6;c W;/ 
-12600c<;;a Rbc-pto-d" + 3840 ui ■ a to -b^^-pc^o -a" 

+ 57&Muj,pUJ,aUJ;bci^:d" - 2104 tJ.pa fee W;/ 
-47520 C^;p Rqab^-.cOJ-.d" - 2520 W;a i?pbgc;d w'^" 

- 1890 W;p Rab;cd + 1890 W;pP i?„fc;cd 

+588a;;pu;;a i?6c;/ + 2352 u;;pa Rbc-.d" 

-5628 cj;p t^:a Rbc ^d + 7728 t^;pa i^tc'^d + 2520 u,p u,a Rb^-.cd 
-8820 W;pa i?fc^;cd + 7920a;;p i?a^6P;cd + 2520 Wjp, Ra\";cd 

-720W;a W;fc W;cW;pd^ - 6480 W ; a ^ ; fc c ^ ; p / + 2 1 60 W ; a W ; W ; c W ; p^d 

+ 3150 Rab:ci^-,p"d + 1800 W;a t^;fcc t^;p^d + 2880 LU -p LU-a LU'^ UJ-bcd 
-3240 U)-a W;p^ W;6cd + 2880 Ul-p Ul-a" LO-bcd + 3600 LO-pa" LO-bcd 
-4500w-p'^a(^;bcd - Q4:80u};pW-a(^;bi^;cd" " 5040 i?pa;6 Wjc/ 
+ 3150 i?ab;p W;cd^ ~ 10080 UJ -a ^ipb to -cd" — 360 

-9000 

^\pab ^\c<f "I" 8730 LO-^abp ^icd" ~l" 9360 UJ-p UJ-^a^;b ^ic^d 

-1260 Rpa;bUJ;c"d + 6300 i?„6;p a;;/d + 19440 uj-a iO;pb io./d 
+5760 UJ-pUJ.ab(^;c"d — 6480 

OJ-^pah d H" 3780 Rpaqb-.c ^■.d^'^ 

+3780 Rpaqb-cUJ-d"" - 1260 Rpaqb;cU)'P''d + 1470 C^;p i?ab;c/ 

+462 u;;p i?afc;/d + 798 u;;p Rab'^cd - 2730 W;p i?a^;6cd 
-1050 W;pi?ga6'';cd'^ -2880 W;aa;;bW;pc/ + 11520 a;;abW;pc/ 
-360 - 1080 + 1080 

-6840 W;afcW;p^cd + 720a;:pPtJ:afccd - 360 CJ;p CJ;a W;6c/ 

+ 10800 3960 

^■,p^;a'^;bc d 

+ 4320 

'^;pa '^;bc d 

+7200a;;pa;;a W;6''cd - 23760 a;;pa Wj^^'cd + lOSOwja Wjcdp*" 
-6840u;;a6 W;cdp'' + 3600 Wja^ipbcd'' + 3600 u;;a W;p6c^d 
-1440a;;a W;p6^cd - 3240a;;a Wjp^bcd + 6120w;p Wjabcd^ 
+2880a;;pa;;a6/d - 360a;;p a;;a6^'cd - 10440 a;;pa;;a^6cd 
+1800w;a Wjbcpd^ - 3240 tj;atj;6cp''d - 3240w;aa;;6cdp^ 

- 720 W;pa6cd^ — 720 LO-pabc' d — 720 IV-pab^cd 
+540 LU. pa'' bed + 900 U;.pP abed - 360 Ul-abped'' 
-360 to -abp/d + 900 LO-abp^cd - 360 Wjabcpd^ 

+900 + 900 

^•,abcdp 

+ 840 OJipOJia^lb^^' Red 

+3150 u;,p Rab'" Red - 3150 uj,p Ra" -.b Red - 840 u:,a uj,b u,pP i?,d 

+420 UJ.pP OJ-ab Red + 1680 LO.p U-a OJ-^b" Red " 4200 LO.pa LO-b" Red 
+1260 UJ-a OJ-pb" Red - 2520 CO-a OJ;p"b Red + 7560 UJ.p UJ-ab" Red 

- 1 1340 UJ.p UJ,a"b Red - 2520 UJ.pab" Red + 630 UJ,pa"b Red 
+ 1890 W;pP„b Red + 1890 LO,abp" Red " 525 UJ,p W'f Rab Red 
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+525 C^;pf Rab Red - 980 UJ,p UJ,a Rb" Red - 770 U,pa Rb" Red 
+ 1680 LU.pUJ. a W-b to- cRd^ + 4200a;;a'^;6^:pc-Rd'' 
+6720uJ,puJ.,ai^:bc Rd'' 

— 10920 UJ,pa OJ;be Rd"" - 10080 UJ,a UJ,pbe Rd" " 5040 U;,p UJ,abe Rd" 
+2520W;aW:6cp i?/ + 3780 W;pa6c i?/ - 1260 W;afccp i?/ 
-1120 W;p W;a i?hc i?/ + 1820 UJ,pa Rbc Rd" - 1764 Rqab^.e Rd" 
-480 LU.p UJ.g Rrba" -R/d'^ - 6480 UJ..P UJ,q UJ-a OJ-b Re'^d^ 
+7Uu,p Rqa;b Re'^d'' + 5250 W;p Rab,q Re" d"" 

-1680 LU, a -.b ypq Rc'^d^ ~ 15360 UJ -.p LU ,a ^ -.qb Re'' d^ 
+34800 u;,pa LO,qb i?c'/ - 240 LO,p LO,q LO,ab Re^d" 
-1680 U),pq UJ,ab Re^d" + 5040 iO,a i^;pqb Re" d" 
-10080 LJ.aUJ.pbq Re" d'' + 20160 UJ,pUJ.qab Re" d^ 
+32040 U,p U,abq Re"d'' - 2520 U;,pqab Re"d'' 
+5040 a;;pa69 i?c^/ + 2520 Wjabp, Re"d'' 
+2016 UJ,p U,a Rqb Re"d'' + 3024 U;.,pa Rqb Re"d'' 
+ 7140 U;;pW;, Rab Re"d'' + 2100 Rab Re"d'' 

-3150 W;p Rqab^'-.r Re"d'' - 16320 LJ.p UJ,q Rrba" Re'' d"" 
-6720 U.pq Rrba" Re^'d'' " 2016 W;p Rgra^.b Re'' d" 
+420 Wjp Rqarb'" Re^'d" - 3990 LO-p Rq^ ra;b Re'' d" 

~3584u;..pu;..a Rqrb" Rc'' d" + 5824c.;;p„ Rqrb" RJ d" 

-420 UJ,p W'P Rqarb RJd" + 420 LO,pP Rqarb Re' d" 
+2688 U-p U}-a Rqbr' Re^ d" ~ 4368 UJ-pa Rqbr^ Re^ d" 
+3584 W;p UJ,a Rq^'rb Re''d" - 5824u.,pa Rq^rb Re^ d" 

— 3276 LO-^p Rqarb;e Rd^"'^ — 3318 Ul-^p Rqarb;e Rd"^^ 
-2912 UJ,p iO,a Rqbre Rd"'"' + 3472 UJ,pa Rqbre Rd""'' 
-672 U.p Rqarb;e Rd''"") /907200 

+gab {24:0cO-pLO-qUJ.ci^;di^''^ ^^''" + 1050 U-p OJ-q U)''^ Red' ^ 

-1050u.pU>.qto'P Red'" - 1050w;pW;,a;'«i?/;d 

+1470 Ul'.p Lo'.q LO'P Re";d - 1440 Uj\p U;\qU!.e Rd^'" 
+480 a;;p UJ.,q LU.c R^" -'d + 225 LJ,p UJ,qUJ.r Re^d"''' 
+ 1995uj-^pUJ-qUJ-r Re' d^'" + 80 Wjp W;c W;d W'^ Wj^'^ 
+1050a;;p Red'^ i^;q" - 1050uj.p ReP'-dU-q" 

— 280 UJ-^cU-dU-p^ UJ-q" — 240 LO-^p OJ-^q OJ'^ U)'" U)-^ed 
+280 UJ.p UJ'P iV.q'^ LJ,cd + 140 Ul.pP Ul.q'' UJ;cd 
+1680oJ-^pUJ-qU)-cLO'" OJ-d^ — 640 LU ■ p Ll) ■ c LO : q" U)-d^ 

— 1800 LO-pc UJ-q'^ LO-d^ — 760 LO-pLO-qLO-c" U-d^ — 3360 LO-pq U-c" UJ-.d^ 
-720 UU.p UJ.q UJ-e W'P OJ-d" + 1620 'uj.,p Rqe^ UO,d" - 240 U.p Rq^-eU-d" 
+360 oJ-pRe^-qW-d" -320uj,pljj''Puj.qeUJ.d" + 1120 UJ -J UJ-qeOJ-d"' 

— 2936 u-p (jj-c^-q^ U-d" + 1280 LO.pe to ■ q^ LO-d" 

+640a;;pa';q oj-.eUj-.dOJ'^" 

-1260 lj,pRqe;dU'P'' + 1050 tJ;p Rcd.qU'P'i - 200uj,eU.,dUJ.pqU;'P'^ 
+3576 Wjp LO-etO-^qd U'^" + 80 UJ-p UJ-^q UJ-^ed ^0'^" + 160 Wjpg U)-ed tJ'^" 
+420u.,pRqerJ''^^''"'' - 900 Wjp W'"'' - 60 W;p'i?/rc;d W'"'' 

+630W;p W;, Red'P" + 360 UJ.,pUJ.,q Red'"^ - 924uj.pUJ.q Re^-d" 
+384 U,p U,q ReP'"d - 360 U,p U,q Re";d'' + 270 W,p W,q RP<ed 
+360 u-p U-e tO'P U-qd" + 420 U-e tO-pP U-qd" — 600 u-p U-e^ U-qd" 
-240 LU,pLU,cU-,dU.,qP'^ + 480 LU ,c U -.pd U ,qP'' + 480 U .,p U .,cd U -.q^"^ 
-240u,pLU,cU'P U.q'ld - S,4:0u,cU-pP LU,q'^d " 720 U;p U;/ U-q"^ d 
+525 LO-^pPcU-^q'^d — 360u-pU-cU-dU-q'^P — 1080U-etiJ-^pdU-^q'^P 
-60 U;p U-ed tO-,q"P + 420 U-ped UJ-.q"^ + 1260 U;p U-q U''^ U-'ed^ 

+3720a;;pW;g«W;c/ - 2340w;pW;q +3360a;;pu;;/a;;c/ 

+750 U,pPqU.,ed" - 3780 U., p U., q U''^ U.Jd - 4140u.,pU,q'' U.,,:Pd 
+4380 u-p U-^q U'P U.,e"d- 5160 W;p W;/ U.,e"d - 420 U-pPg U.^e"d 
+480 u-p U-qU-etO-d"" - 4320 U-pU-qeU-d^'' + 1500 U-peqU-d"" 
-2040 U,pU,qcU,d"P - 1260u,peqU,d"P '- 240 U ,p U ,q U ,c U'^'^ d 
-360 U-^etJ;pqU''P"d + 2280 U-pU.qeU'P"d + 210 U-pqe U'"" d 
-480 u-p U'^U-qed" - 840 U-p^U-^qed" - 60 U-p U'^ U-ge" d 
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+240 c^;p, cj^q/P - 480 cj.,p uj^cC^-./d" + 240 u;.,p, c^;// 

+360 LU,p LU,c UJ.qP'^d - 180 UO.,p, UJ.,gP'id + 180 UJ..P LU-'P UJ.,g\d 

+630w;pf a;;,«c<i - 840w;p W;cW;g«/ + 960a;;pcW;g«/ 
-240w;pa;;cW;q«Pd + 360 W;pc + 180w;p a;=Pa;;cdq« 

+630 a;;/W;cdq'^ +2220 +480 W;p5W;c/'^ 

-4440 W;p W;, W;// - 360 a;;p, wi/d" - 480 W;, u;;/«d 
+120 W;p, a;;/«d - 480 W;p t^;c + 660 W;^^ ^Jq'^ 

+ 1740 W;p cj., cj^P«,d + 240 uj,pg oj--P\d - 360 u^,p uo-.q^J"" 

-mQOJ.pOJ.qcd"" - 120 W;pW;,cV - 120a;;pW;gc^'d 

+240 uj,p a;;ge^/ + 60 io,p uJ-^qc^^a + 240 io,p uj.gP^d" 
+240 u,p uj,qPc''d - 180 W;p uj,qP\d + 540 a;;p 

+240W;pW;,«C^d - miO,pUJ,q'iPcd - 180W;pW;cd/« 

+540 uj.p uj-odq'^P + 360 LO.p LO.cd^q'i + 240 a;;p W;/,/ 

-180W;pW;c%'^d - 180W;pW;c^d,« " 70 W ; p W ; 5 W ^ " W ^ ' iJ^d 
-140 W;p Red + 175 W;pP W;,^ R^d " 560 W;p LO'^'^ Red 

+ 140 W;p5 u!'P'''Rcd + 420 a;;p LO.q''P Red - 700 W;p W;g a;;c i?/ 
-700 u)\p iv-c u};q'^ Rd^ + 350 uj-pc i^J-q'^ Rd^ + 420 a;;p W;, w'^ -Rd'^ 

-\-7QQuJ.pJP U-qeRd"^ - 560ui.pU.ciO-qP Rd^ +'i80ui.pU).qU}.c^ Rd"^ 
~UQoJ..pqOJ..P' Rd'^ ^MQ0J..pij..qP Rd'' - ^'20uJ:pOJ:qPcRd'' ' 
-420 LU-p OJ.Pq Rd'^ - 40 UJ-p LO-qLO-cOJ-d R^'^ " 1640 LO.p UJ-cOJ-qd R^'^ 

+920 w,pc ^.qd - 240 w,p u,q W;ed ii^" + 420 uj,p u,qcd R''" 

+540 IV-.p LO-cdq i?^' + 210 LO.p UJ.q Red R^'' " 840 LO-p UJ-q/ RJ'd'^ 

+ 1800 a;;pW;g«^i?e^/ + 3820 w.pw.qw./i?^'/ 
+700a;;pqW;/i?c«/ + 840w;pW;,/ i?>/ 
+2680 w.p w.q w./ i?c'^d'^ + 920 io.pq w./ i?^'^^'^ 
+MQuj.pUJ.qPr Rc'^d' + 864a;;pa;;, i?> R^'^J' 
+700uj'.pU)]qUi.,rU>''' R/d^ - IQs'uj.p UJ -q Rr'^ R/ d^ 

-2100 U>..p UJ.q OJ-r tO'P Rc'^d'' - 280u).pU}'Pu).qrR/d'' 

+960w.pW.qi?/«/i?// - 1560u).pU).qRr''s^ Rc^d'' 
-3280 cL;p i?/"'^ - 280 uj..pUj..qrc Rd^'^'' 

+2480a;;pW;,cr-R/'^'^ - 1176w;pW;g i?„i?/«'^ 
+324 a;;p W;, Rrsc" Rd"""' + 1120 w.q^r Rd'""' 
+432 w.p w.g Rrc Rd'^^'' - 320 w.p w-c w.,^ 

+400 w'pc U.qr Rd''^'^ - 1360 W;p UJ.^cr Rd''"'' 
-80 Wjp W;, Rrsc" Rd"" + 40 W;p W,q Rr\c Rd''"' 

-792 Lu,pLu,qRrs,J^ Rd'P"- 

+560 u;;p uj.^q Rr^'sc Rd'"l /302400 

+9ab9cd (-12 W;p W;, a;;r uj''P w''' w''' 

+108 a;;pW;ga;;,. -60 u;;pU;;gW'Pa;'«a;;/ 

+ 21 Ul.^p UJ'P UJ .^q'^ UJ +35 W;pP LO .^f^ LO — 12 LO-p LO-q LO-r' ^'""^ 
+84:CO.pqU}.r' ^'""^ — 588c0.pCO.qCO.r<^''' <^'"'' — 3QQ iO-p LV-q LV-r'^ Ul'P^ 
+64:U)-pqU}-r'' '^'"'^ + 252 U)-p U)-q U)-r 1^'" 1^''''' — 12 UJ -p UJ'P UJ -qr UJ''^'^ 

— 36 UJ.p UJ.q UJ'P UJ-r'''' + 252 UJ -p UJ ■qUJ''' UJ + 252 UJ -p UJ ■q'' UJ ./P 

— 216 UJ-p UJ-qUJ'P UJ-/'^ + 2>06uJ-pUJ-qP UJ-/'^ — 108 U -p UJ -q UJ -r ^jJ'"'''' 
+ 108w.pW.qrUj'^P''''> + 108uJ,pW.qrUj'^'''''P - 36 W;p W;, W;/"^ 

-18 UJ,pUJ,qUJ,r"'''' + 162 UJ.,pUJ.,qUJ,/"'' + 5A UJ ,p U ■qUJ'P'^ / 

+ 126 w.p UJ.q uj./' RP'^ + 126 ^.p w., i^P*^ 

-108W;pW;gW;^W'Pi?«'' 

+306 w.p w.g w./ i?**^ + 36 w.p w.g uj-rs R"'"") /181440 (Ble) 
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APPENDIX C: EXPANSION TENSORS FOR Gfl„ 



We present the explicit expressions for the expansion tensors for Gfin- 



k2 



g(°) = _:_ (Cla) 

tin g g2 til ^ ' 

'^fln.ah = 72 e^"" {Sw-aW-b - 'iw-ab -2w-pW'P Qab + W-/ gab + Rab) 

^fin abc = 1/1/1 2m) (8W:a W;fcW;c - flafcic - 12W;a W;bc + 2W;a{,e - 4W;pW;a W'^gbc 



144 

2 W;a W;p^ + 4 W;a'' Qbc - W-p^a 9bc + 2 W-a Rbc) 



180 e4,„ V- . aoSl+2e'- W..a SI 51 - W,a9bc) (Cld) 

^fin,a6c = ^^^^^-^^ (192 W;aW;bW;cW;d - 60 W;ai?fcc;<i - 576 W;aW;6W;cd 

+ 144w.<jb W-crf + 18i?ah;cd + '^^^.W-aW-bcd " ^Aw-abcd 

+ 80 W;a W-bRcd - 40 Wjab -Red + 5 i?ab -Red + 4 i?pa9fc -Re^d"^ 

-25q6 {72 W;pW;cW;dW'P - 15W;pi?ed'^ + 15W;pi?e^;d 

— 36 W;e W-dW-p^ — 36 W;p W''' W-cd + 18 W;cd 
-144?«;p W;cW;d^ + 36W;pcW;/ " GW-cW-pJ" 

+42 w-c w-p^d - 18 w-p w-cd^ + 54 W;p w-/d 

+ 12W;pe/ - 3W;pe^d - 9W;p^cd 

-9 w-cdp^ + 10 w-p w'P Red - 5 w-pP Red 

+ 10 ?Z;;p W;e i?d^ - 5 W;pc -Rd^ " 10 W;p W ; q Rc^d'^ 

-10u;;p,i?eV} 

+.9a65cd {12 W;pW;gW'Pw'' - 14W;pW'PW;q' + 5w.pPw-q'^ 



-28w.pW.qW''P'^ +Aw.pqW''P'i 

+12 m;;p w,g''P + 6 u;;p W;, 



4320 ^^'^^ ^'^'^ + 12 «^;c - 4 F^^^e 

+ 8W;aW;6(5e <5d -4:W.abSl Sd+^l Rab} 

-gab {28 w-cW-d- 8w.cd + Rc 



4e2- (2w;pw;'f-w;;pf) <S^<5Si?ed} 



+ gabgcd {3 W;p - W;p^}) 

Y^^^^ (16 - 51 SI SI - 12 - 5^ gab + 5ab 5cd) (Cle) 
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